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We present a calculation of the properties of vibrational states in deformed, axially-symmetric 
even-even nuclei, within the framework of a fully self-consistent Quasparticle Random Phase Ap- 
proximation (QRPA). The same Skyrme energy density and density-dependent pairing functionals 
are used to calculate the mean field and the residual interaction in the particle-hole and particle- 
particle channels. We have tested our software in the case of spherical nuclei against fully self 
consistent calculations published in the literature, finding excellent agreement. We investigate the 
consequences of neglecting the spin-orbit and Coulomb residual interactions in QRPA. Furthermore 
we discuss the improvement obtained in the QRPA result associated with the removal of spurious 
modes. Isoscalar and isovector responses in the deformed 24_26 Mg, 34 Mg isotopes are presented and 
compared to experimental findings. 



I. INTRODUCTION 

The response of many-body systems to an external, 
weakly coupled field provides much insight regarding the 
correlations existing among the particles composing the 
system, and the forces acting among them. 

While the most familiar application of mean field the- 
ory is to describe stationary states, its extension to time- 
dependent states provides the basis for a theory of small 
amplitude oscillations known as Random Phase Approx- 
imation (RPA) for normal systems and QRPA for super- 
fluid (superconducting) systems displaying Quasiparticle 
excitations. 

Especially, small amplitude oscillations describe the 
dynamics of the nuclear surface. Couplings to the surface 
will influence the quasiparticle motion, renormalizing the 
effective mass (w-mass m u ), in turn leading to couplings 
among the different vibrational modes. 

A textbook example of such couplings is provided by 
the breaking of the giant dipole resonance in deformed 
nuclei (inhomogeneous damping [H-0), in keeping with 
the fact that a permanent deformation can be viewed as 
a quadrupole vibration of finite inertia and of vanishing 
restoring force. Within this context we refer to Fig. 10 
of ref. ,4j , in particular to the two-peak photoabsorption 
cross section (a(E;El)) of the 150 Nd nucleus, marking 
the onset of static deformation in the Neodynium iso- 
topes as a function of mass number. 

As the smooth increase of the FWHM of cr(E; El) as a 
function of A indicates, this cannot be a yes-or-no effect. 
In fact, the FWHM of <r(E; El) associated with the tran- 
sitional nucleus 148 Nd is not very different from that of 
150 Nd. Within this scenario, quadrupole deformations, 



static or dynamic, will also modify not only the proper- 
ties (centroid and width) of the GDR but also those of 
the GQR and, in deformed nuclei, that of the GMR. The 
same line of reasoning also applies to other multipolari- 
ties of the static and dynamic deformations of the mean 
field, which should not necessarily be only quadrupolar. 

Because exotic nuclei, in particular neutron halo nuclei 
are, as a rule, more polarizable than nuclei lying along the 
stability valley, one expects this to affect the modes and 
the associated renormalization to be especially important 
for such nuclei. Here, one may mention the pygmy reso- 
nance [f| and the inversion of the usual sequence of single 
particle energies [6] . A consistent treatment of the vibra- 
tional modes is the first step on the way to address such 
properties. 

The low-frequency collective excitations are quite sen- 
sitive to the shell structure near the Fermi level as well 
as to the nuclear surface shape, and one expects that 
new kinds of collective excitation will emerge under new 
situations of nuclear structure. To investigate such pos- 
sibilities, many calculations have been made using the 
self-consistent RPA based on the Skyrme-Hartree-Fock 
(SHF) method 0,1] and the Quasiparticle-RPA (QRPA) 
including pairing correlations 043 • A number of simi- 
lar approaches using different mean fields have also been 
carried out [Ia - l22l ]. 

Recently new iterative methods have been developed to 
calculate RPA strength functions for both spherical (23[ 
and deformed [24l l25j nuclear systems. Low-frequency 
RPA modes in deformed nuclei close to the neutron drip 
line have been studied |2(| [53] , taking also into account 
pairing correlations [HI, [2t| ■ These latter calculations are 
based on a BCS approximation which does not take into 
account continuum coupling effects, typical of drip-line 
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nuclei. A proper theoretical description of such weakly 
bound systems requires a careful treatment of the asymp- 
totic part of the nucleonic density. An appropriate frame- 
work for these calculations is provided by the Hartree- 
Fock-Bogoliubov (HFB) formalism, solved in coordinate 
representation [30|, HH especially for spherical nuclei [HJ , 
or more conveniently in the configuration-space approach 
for deformed nuclei (33j . 

A quantitative description of excitations in exotic nu- 
clei is given by fully consistent QRPA calculations on 
top of an HFB ground state, such that the same effective 
interaction is used for both calculations. A fully con- 
sistent HFB+QRPA approach with the Gogny effective 
interaction for spherical and deformed nuclei has been 
developed in a harmonic oscillator basis 34]. Standard 
QRPA equations have also been solved in a cylindrical 
box with the Skyrme effective interaction, not including 
neither spin-orbit effects nor the Coulomb residual inter- 
actions [35J, l36| . 

Within this context we discuss in the present paper 
a consistent approach to describe linear response in de- 
formed nuclei within the framework of HFB+QRPA. Sec- 
tion HD discusses the elements used to work out a software 
to implement such a program. In Sec. IIII1 we provide 
detailed information regarding Hartree-Fock-Bogoliubov 
ground states. Section Hill also illustrates some basic re- 
sults of the QRPA software developed by us to treat de- 
formed nuclei when applied, for a consistency check, to 
the case of a spherical system. In this Section we also 
carry out a discussion concerning the spurious modes. In 
Sec. |IV]the response functions of 20 O, 24 ~ 26 Mg and 34 Mg 
are shown and discussed in comparison to available data 
and to other calculations. Conclusions are drawn in Sec. 



II. METHOD 

As the first step in the self-consistent calculation of 
excitations in axially deformed and reflection symmet- 
ric nuclei, HFB equations are solved in a finite harmonic 
oscillator (HO) as well as in a transformed harmonic os- 
cillator (THO) basis. In both cases, a discretization of 
the positive energy continuum is carried out. We use the 
new version (known as (101)) of the program HFBTHO 
[Hj], choosing for 20 O, 24 ~ 26 Mg and 34 Mg, a quasiparti- 
cle energy cutoff E cut — 50 MeV, and set N s h = 15 HO 
(THO) shells. This code allows to perform HFB calcula- 
tions with arbitrary Skyrme functionals together with a 
density-dependent pairing delta interaction ]39J 
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Poo ( r ) being the associated isoscalar density and P a the 
spin exchange operator. In the following, we use the 
Skyrme functional SkM* [40| and a pairing interaction 
with the same parametrization as that adopted in ref. 
Fill . i.e. for 20 O and 24 ~ 26 Mg we adopt the parameters 



V = -280 MeV fm 3 , Vi = -18.75%, 7 = 1 while for 
34 Mg Vo = -295 MeV fm 3 , V x = - 18.75 V , 7 = 1 are 
employed. It corresponds to a mixed surface- volume type 
of pairing potential. 

We diagonalize the HFB Hamiltonian in configuration 
space. We then diagonalize the density matrix p ob- 
taining the canonical basis. We have checked the nu- 
merical accuracy of the procedure comparing the quasi- 
particle energies obtained in the original diagonalization 
with those obtained diagonalizing the HFB Hamiltonian 
in the canonical basis. We have found that the values 
agree within 10 -4 MeV. Finally, we switch to coordinate 
space, tabulating the canonical wave-functions together 
with their first and second derivatives in a grid of 30 x 
30 (Gauss-Hermite) x (Gauss-Laguerre) points for Gauss 
quadrature integration in cylindrical coordinates (r±, z). 

The QRPA basis is obtained using pairs of these 
canonical wave functions such that the dimension of 
the QRPA-Hamiltonian matrix does not exceed the size 
20000 x 20000. For this purpose, one first omits the 
canonical states i that have single-particle energies Ei 
greater than some £ C rit- A second cut is made exclud- 
ing those QRPA quasiparticle pairs displaying occupa- 
tion probabilities less than some small v^ rit or larger than 
1 — v 2 crit . In the following calculations these parameters 
are given the values e cr n = 200 MeV and v cr i t — 10~ 2 . 
The QRPA excited states |A) are described in terms of 
the quasi-boson operators 



E 
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^KK' a K a K' ~ YKK' a K' a K 
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acting on the QRPA correlated vacuum |0) (Qx\0) = 0). 

The operators or K , olk are the canonical quasiparti- 
cle creation and annihilation operators respectively and 
Xkk' > ^kk> are t ne amplitudes of the two quasiparticles 
excitations {K, K'}. The matrix elements between differ- 
ent QRPA basis states {K, K'} and {L, L'} are expanded 
respecting the selection rules of the vibration's quantum 
numbers f2 and n. That is, f2 = £Ik + Qr 1 = + 
and 7r = ttk ' ^K' = 77 l 1 n L'- Here, £1 is the projection of 
the angular momentum on the symmetry axis z and n is 
the parity. Details concerning the calculation of matrix 
elements are given the Appendix. 



III. BASIC RESULTS 

A. HFB ground states 

First, potential energy curves are calculated as a func- 
tion of the deformation parameter /3, defined as 



n <Q> n + <Q > p 



(3) 



The quantity < Q > q is the average value of the 
quadrupole-moment operator Q = 2z 2 — r\ for protons 
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(q = p) and neutrons (q = n). The QRPA calculation will 
be performed based on the HFB solution corresponding 
to the absolute minimum of the potential energy. 

Fig. Q](left panel) shows the potential energy curves for 
the nucleus 24 Mg, comparing the results obtained in the 
HO and THO basis. One finds a pronounced minimum 
corresponding to the prolate deformation ft — 0.39. Sim- 
ilar HFB calculations have been carried out for the other 
isotopes discussed in the following. The nucleus 20 O is 
found to be spherical, while 26 Mg is oblate (/3 = —0.18) 
and 34 Mg is prolate ((3 = 0.36). The ground state prop- 
erties of the four mentioned nuclei in the HO and THO 
basis are summarized in Tabic [I] 

The right hand panel of Fig. [T]shows for 24 Mg the pair- 



ing energy E, 



pair 



^Tr(AK) calculated separately for 
protons and neutrons, where k is the expectation value 
((HFB\P+\HFB) = (HFB\P\HFB)) in the super- 
fluid ground state (\HFB)) of the pair addition/removal 
(P+ = Y.v c t c t I p = X^cpc^ c t'- single-particle 
creation operator) operator, that is of the pair field (ab- 
normal density) . A stands for the functional derivative of 
the energy E\p, k] with respect to the abnormal density 
(pairing gap) [H,|4l|. 

Returning now to Fig. [TJ one may notice the differ- 
ence between the curves for protons and neutrons. The 
Coulomb field influences the density of the nucleus and 
so the pairing interaction. At /3 = 0.39 the total pair- 
ing energy is zero, so QRPA calculations reduce to RPA 
ones where the excitations are only in the particle-hole 
(ph) channel. Full QRPA calculations are performed in 
the other three nuclear systems 20 O, 34 Mg and 26 Mg for 
which the HFB ground state displays neutron and proton 
pairing correlations respectively. 
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FIG. 1. HFB potential energy curves Ehfb (left hand panel); 
neutron (thick lines) and proton (thin lines) pairing energies 
E pa i r (right hand panel) . The quantities are calculated in HO 
(dashed curves) and THO (solid curves) basis as functions of 
axial deformation parameter /3 in the nucleus 24 Mg. 



B. Response for the spherical nucleus 48 Ca 

To test our deformed QRPA code developed on the 
basis of the deformed HFBTHO program by M. Stoitsov 
331 . we compare our results to those obtained with the 



spherical QRPA code developed by J. Terasaki et al. (JT) 
14], making use of the same Skyrme interaction. Fig. 
[2] shows, for the spherical nuclear system 48 Ca, such a 
comparison with calculations based on either (i) spheri- 
cal basis states with a hard-wall boundary condition at 
20 fm, and (ii) the two types of deformed basis, HO and 
THO. The response functions of Fig. [2] have been calcu- 
lated with the help of Eqs. (1), (2) of ref. [42|. In the 
following, we use strength functions defined as 



r/2 



* (E-E x y+ry4: 
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for the multipole operator Fj n . If not specified in the 
text, the value F = 1 MeV is used as in ref. [35] in the 
calculation of the IS monopole and quadrupole transition 
operators 
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and for the IV dipole operator 
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The agreement between the strength functions labeled JT 
42] with those labeled HO and THO shown on Fig. [His 
actually quite rewarding, in keeping with the completely 
different routes taken in the calculations. In all the three 
cases the strength function exhausts about 98 % of the 
Energy Weighted Sum Rule (EWSR). 

For the 2 + mode in 48 Ca, the degeneracy between the 
components fi" = + ,±l + ,±2 + is reached up to the 
order of 10" 3 in both the HO and THO basis. 



C. Self-consistency and spurious states 

Let us now turn to the spurious solutions of the QRPA 
equations in the dipole modes with Q n = 0~, ±1~ and in 
the quadrupole modes with fl* = ±1 + . Spurious modes 
are, in the present case, due to the translational and 
rotational symmetry breaking, respectively, of the HFB 
ground state [4l|. The quadrupole modes with fi^ = + 
also contains spurious states associated with particle- 
number nonconservation. Former works on QRPA both 
for spherical nuclei [3] and for deformed ones [34[ show 
that the spurious states which should be at zero energy 
[43| have instead values from a few keV to 1-1.8 MeV, 
well separated from, in any case, (physical) states. These 
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TABLE I. Ground state properties of 20 O, 24 ' 26 > 34 Mg obtained by the deformed HFB calculation in the HO and THO basis. Chemical 
potentials A 9 , deformations /3 q , average pairing gaps A 9 , root-mean— square radii •/ '{r 2 ) q for neutrons (q = n) and protons (q = p), and 
the total binding energies Ejjpg are listed. 
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FIG. 2. Isoscalar (IS) and isovector (IV) 2 + strength functions 
in 48 Ca in the spherical basis of J. Terasaki (JT, dashed-dotted 
curve), in HO (dashed curve) and THO (solid curve) basis. 



nonzero energies may be taken as a measure of numer- 
ical accuracy with which linear response is calculated, 
and are strongly related to the choice of basis states and 
truncation of basis size. 

The spurious states obtained in our fully consistent 
calculations lie either on the real on the imaginary axis, 
and their absolute value never exceeds 700 keV. In Fig. 
[3] and 2] we show results for 24 Mg, comparing the self- 
consistent results with those obtained (i) neglecting both 
Coulomb and spin-orbit residual interaction and (ii) ne- 
glecting only the Coulomb residual interaction. It is seen 
that the absolute values of the energies are progressively 
reduced going from case (i) (about 3 MeV) to case (ii) (2- 
2.2 MeV) and finally to the self-consistent solution (less 
than 0.5 MeV). One also notes that in case (i) the en- 
ergies of the spurious states lie on the real axis, while 
in case (ii) they are found on the imaginary axis. The 
fully consistent QRPA calculations give imaginary ener- 
gies for the dipole spurious modes and real energies for 
the quadrupole spurious modes, all below 500 keV. It is 



satisfactory that the energies of the spurious modes get 
closer to zero as more terms of the residual interaction 
are included. We ascribe the remaining distance to zero 
to the truncation of the basis and to numerical inaccura- 
cies. Within this scenario one may posit that a spurious 
mode obtained by fully consistent QRPA calculations at 
"imaginary energy" 100-500 keV is as good a solution 
as that corresponding to a spurious real energy mode at 
100-500 keV. 

A similar study has been performed on the energies of 
the dipole spurious modes of 20 O. Also for this nucleus 
if both Coulomb and spin-orbit parts are neglected, the 
spurious energies are at 2.5-3 MeV. The values jump to 
about 1.5-2 MeV on the imaginary axis if one adds the 
spin-orbit contribution to the residual interaction. Fi- 
nally, the energy of a spurious state obtained with a fully 
consistent QRPA calculation can be real or imaginary 
not exceeding 600 keV. For the deformed nuclei 26 Mg 
and 34 Mg we obtain from the fully consistent QRPA cal- 
culations with SI* = 0~~,±1 - and fl n = ±1+ spurious 
states with either real or imaginary energy, in any case 
the modulus not exceeding 700 keV. For 20 O and 26 ' 34 Mg 
the quadrupole spurious modes with Sl n = + are all real 
at around 1.8-2 MeV. We ascribe this problem concern- 
ing the SI* = + modes mainly to the rather crude cut 
in the occupation probability given by v cr it — 10~ 2 . We 
have checked in the case of 26 Mg that reducing the value 
of v cr it to 10~ 3 brings the energy of the + spurious state 
from 1.8 MeV (on the real axis) down to 1.1 MeV (on the 
imaginary axis) . The transition strength is only affected 
for energies below about 3 MeV. 

Similar investigations of the spurious modes for calcu- 
lations with the Gogny force have been reported by S. 
Peru et al. in references [44[ and [34[ for spherical and 
deformed shapes, respectively. For spherical nuclei, it is 
found that the spurious modes lie at very low energy, 
about 3-5 keV, when all terms of the interaction are in- 
cluded, whereas they may move up to about 2 MeV when 
leaving out parts of the interaction. For deformed nu- 
clei [3J|, the result is qualitatively the same as in the 
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present work. However, the Gogny results are better 
for the = + mode, which comes at a very low en- 
ergy, but worse for the translational invariance modes 
Cl n = — ,±1 — , where the spurious mode energies can 
come as high as 1.8 MeV. 
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FIG. 5. IS 2+ strength lunctions in 20 O, in HO (upper panel) 
and THO (lower panel) basis without the residual spin-orbit (SO) 
and Coulomb (C) interaction (dashed curve) and with all the terms 
included (solid curve). 
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FIG. 3. Energies on the complex plane of the spurious state for the 
1 — mode with Q n =Q~ (circles) and Q n =l~ (crosses) projections 
in 24 Mg in HO basis. The figure illustrates the change in energy of 
spurious modes when specific terms of the interaction are omitted. 
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FIG. 4. The same as Fig. [3] but for the 2+ mode with fr=l+ 
projection. 



IV. STRENGTH FUNCTIONS 
A. 20 O 

The IS quadrupole 2 + mode is shown in Fig. \5\ cal- 
culated with the HO as well as with the THO basis 
states. A large fraction of the strength function is concen- 
trated in three energy regions: a low-lying mode around 
3 MeV, a small peak around 10 MeV, and a giant res- 
onance around 20 MeV. The two calculations display a 
remarkable agreement with respect to the strength func- 
tion below 15 MeV. Whereas the calculations agree on 
the centroid and the strength of the giant resonance, they 



display differences with respect to the splitting of the two 
components of the resonance, and on the distribution of 
strength between the two peaks. When we omit both the 
spin-orbit and Coulomb residual interaction, in the giant 
resonance region the strength functions are shifted down 
in energy about 0.3 MeV for both the basis used, and the 
peak heights increase by 15% for HO and 25% for THO. 
Concerning the two peak energies below 10 MeV, they 
are both shifted up by about 0.2 MeV while the peak 
heights are lowered by 20% for both the two basis. 

Fig. [6] shows the IV giant dipole resonance in the en- 
ergy range (15-25) MeV. In the THO basis the peak is 
at 20.2 MeV, while that in the HO basis is shifted down 
about 1 MeV. 

The strength functions calculated without the spin- 
orbit and Coulomb interactions can be compared to those 
obtained by K. Yoshida et al. [35|, who did not include 
these terms. We find a very good overall agreement. On 
closer inspection, the energy of the giant resonance peaks 
in our calculation in the THO basis is about 0.3 MeV 
higher than in ref. [35]. 
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FIG. 6. The same as Fig. [5] but for IV 1 strength functions. 
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B. 24 < 26 Mg 
B.l. IS 2 + modes 

In Figs. [5] we display the response functions of the 
IS quadrupole mode of the deformed prolate and oblate 
24,26 Mg isotopes. 

For both these nuclei, the first low lying state is 
found to have projection — 2 + , corresponding to a 
7-vibration, while the /3-vibration with Q — + is sit- 
uated at a considerably higher energy. Comparing the 
nuclei, one sees that 24 Mg, which may be described as 
a deformed closed shell nucleus for the prolate shape, is 
more stiff towards vibrations than the oblate 26 Mg. This 
is all in qualitative accordance with data [451 ]. However, 
a more careful comparison of calculations and data reveal 
some discrepancies, especially for 26 Mg. 

In 24 Mg, the calculated energy of the lowest = 2+ is 
4.14 MeV, to be compared to the experimental value 4.23 
MeV of the gamma-band head, the state 2^" (45|. Like- 
wise, the calculated energy of the lowest i7 x = + is 7.08 
MeV, should be compared to the experimental value 6.43 
MeV. Both these states are not very collective, composed 
equally of proton and neutron particle-hole excitations: 

for the 
for the 

= + state, respectively. Still, the RPA root is shifted 
down in energy by about 1.5 MeV relative to the single 
ph excitations. These assignments are in reasonable ac- 
cordance with experimental information from 7— decay 
[45l | . particle transfer [if3], (e,e') inelastic scattering (47j . 
(a, a') scattering Q and (n,Tr') scattering [49 ], 

In 26 Mg, the lowest J7 7r = 2 + excitation is calculated 
to be at 1.31 MeV, considerably below the experimental 
state 2^ at 3.0 MeV, while the lowest 0,* = 0+ exci- 
tation is at 2.75 MeV, to be compared to the state 0^ 
at 3.64 MeV. Compared to 24 Mg, the fT = 2+ excita- 
tion is predicted to be more collective. This is not in 
accordance with the experimental B(E2) values from the 
ground state, which are of about equal magnitude for the 
two nuclei. Also, the considerable matrix element for ex- 
citing the 17*' = + state predicted by the calculations is 
in disagreement with the long life-time seen experimen- 
tally of the O2 state. 

The rather low collectivity of the j3- and 7-vibrations 
follows a general trend. With a permanent deformation, 
most of the quadrupole collectivity is tied up in the de- 
formation of the mean field. This will in turn influence 
the renormalization effects these modes will have on the 
single particle motion. In this context, one can quote the 
rather different pattern observed in the distribution of 
matrix elements contributing to the induced pairing in- 
teraction in spherical and deformed nuclei, respectively 

0, [Hi- 
Next, we turn our attention to the giant vibrations sit- 
uated in the energy region 15-25 MeV for all cases, with a 
characteristic splitting between the different projections 
Q. For the prolate nucleus 24 Mg, the SI 71 ' = + vibration 
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TABLE II. Theoretical mean energy values (in MeV) obtained 
with the Gogny force by S. Peru et al., with the SkM* and 
SLy4 Skyrme force in HO and THO basis (present work), and 
experimental mean energy values by Youngblood et al. and Ir- 
gashev et al. of the IS giant monopole, quadrupole resonances 
(calculated in the energy interval [i? m m, -E max ] = [9, 41] MeV) 
and IV giant dipole resonances (calculated in the energy in- 



terval [Emi n , -Ema: J 



[10, 29] MeV) in 24 Mg. 



is along the longest axis and acquires the lowest energy. 
For the oblate nucleus 26 Mg, the opposite behavior is ob- 
served, and here the splitting is less pronounced due to 
the smaller value of the deformation parameter Q. 
In Tables [Til IHH the mean energies 



E(GR) 



mi [E, 



nun? ^max] 



™0 [-^min, ^max] 



(8) 



with moments 



mx [£min, £ max ] = H E k 1°) 



(9) 



calculated for the 2 + modes of 24 > 26 Mg in the energy 
range [E m i n , E max ] = [9,41] MeV, are given in HO and 
THO basis. First, one sees that the mean value in the 
giant resonance region is practically not affected by the 
choice of basis. Secondly, one may address the 1 MeV 
difference between the results obtained with the SkM* 
Skyrme force and those with the Gogny force, by S. Peru 
et al. [34J. Part of this difference may be caused by the 
lower effective mass m*/m — 0.7 for the Gogny force, as 
compared to m*/m = 0.8 for the SkM* force. For this 
reason we performed for 24 Mg a calculation also with a 
SLy4 force where m* /m — 0.7. The centroid energy then 
increases by roughly 0.6 MeV, indicating that about half 
of the difference between the two calculations can be as- 
cribed to the effective mass. Finally, the results may be 
compared to the experimental mean energy of 16.9 ± 0.6 
obtained by Youngblood et al. J52], and one sees that 
all the theoretical calculations overshoot this value by at 
least 2 MeV. 

The theoretical results obtained in the present work 
with the SkM* and SLy4 forces exhaust 83-84 % of the 
EWSR exceeding the experimental result by about 20 %. 
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26 M E 


IS 


IS 


IV 


E(GR) 


r = o + 


«T = 2+ 


J n = 1~ 


IirtsrvEil 


19 411 


19 411 


no 291 


Theor. D1S (Peru et al.) 


22.0 


21.0 


22.9 


Theor. SkM* (HO) 


22.1 


20.3 


20.3 


Theor. SkM* (THO) 


21.9 


20.3 


20.2 


Exp. (Fultz et al.) 






20.6 



TABLE III. Theoretical mean energy values (in MeV) ob- 
tained with the Gogny force by S. Peru et al., with the 
SkM* Skyrme force in HO and THO basis (present work) 
of the IS giant monopole, quadrupole resonances (calculated 
in the energy interval \E m i n , -Emax] = [9, 41] MeV), and of the 
IV giant dipole resonance (calculated in the energy interval 
{E min , £ max ] = [10, 29] MeV) in 26 Mg. For this last resonance 
the experimental mean energy value by Fultz et al. is also 
given. 
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FIG. 7. IS 2 + strength functions in 24 Mg, in HO (upper panel), 
and THO (lower panel) basis for the Vl w =0 + (solid curve), 1 + 
(dashed curve) and 2 + (dashed-dotted curve) excitations. 



This difference in the EWSR as well as in the mean en- 
ergy between the present work and the experimental data 
of Yovmgblood et al., is apparent from Fig. [5J which 
shows the comparison between calculated and experimen- 
tal strength functions. The shape of the experimental 
curve is reproduced in its overall features, in particular 
for T=3 MeV. However, according to the calculations, the 
central peak around 20 MeV is too pronounced, and there 
is too little strength towards lower energies. It is of no- 
tice the sharp peak experimentally observed at 15 MeV. 
To which extent it may be connected with the £l n = + 
mode which in our calculations appear blue shifted by 
2.5 MeV (i.e. at 17.5 MeV) is an open question. 

Fig. [TU] shows for 24 Mg the effect of leaving out the 
spin-orbit and Coulomb parts of the residual interaction. 
One sees that this would lead to a downward shift of the 
giant resonance by about 0.9 MeV. Comparing to the 
work of K. Yoshida et al. (3f| one finds that their peaks 
are shifted further down by about 0.9 MeV with respect 
to those obtained in the present work without the spin- 
orbit and Coulomb terms. This remaining shift should 
be compared to the equivalent shift of 0.3 MeV discussed 
above for 20 O. It should be ascribed to the differences 
between the two calculations, the renormalization of the 
interaction in ref. [35;], and the different basis used. 



B.2. IV 1' modes 

Figs. [TT] and [12] show the response functions of the 
IV dipole modes, of the deformed prolate and oblate 
24 ' 26 Mg isotopes. The IV giant dipole resonances show 
a two-peaked structure. For both nuclear systems the 
low-lying part of the resonances is given by a defined 
peak at around 16 MeV and 18 MeV for 24 Mg and 26 Mg 
respectively. The higher energy part of the strength is 
fragmented, especially for the THO basis, in fair agree- 
ment with the responses given by S. Peru et al. [34| . In 
the HO approach one may see a two-peaked structure up 
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FIG. 8. IS 2 + strength functions in 26 Mg, in HO (upper panel), 
and THO (lower panel) basis for the Vl w =0 + (solid curve), 1 + 
(dashed curve) and 2 + (dashed-dotted curve) excitations. 
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FIG. 9. Fractions of IS quadrupole EWSR in 24 Mg, in HO (up- 
per panel) and THO (lower panel) basis. Dashed and dot-dashed 
curves are obtained from folding of QRPA spectra with a Lorentian 
distribution having T=l MeV and T=3 MeV respectively. These 
are compared with the experimental data by Youngblood et al. 
(solid curve). 
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FIG. 10. IS 2+ strength functions in 24 Mg, in HO (upper panel) 
and THO (lower panel) basis without the residual spin-orbit (SO) 
and Coulomb (C) interaction (dashed curve) and with all the terms 
included (solid curve). 



FIG. 12. IV 1 — strength functions in 2e Mg, in HO (upper panel), 
and in THO (lower panel) basis for the f2 7r =0 - (solid curve), 1 — 
(dashed curve) excitations. 



to around 26 MeV for 24 Mg and a denned peak at 21-22 
MeV for 26 Mg. 
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FIG. 11. IV 1 — strength functions in 24 Mg, in HO (upper panel), 
and in THO (lower panel) basis for the Q w =0~ (solid curve), 1 — 
(dashed curve) excitations. 

The fraction of the EWSR in 24 Mg for the IV 1" mode 
calculated in the present work for two values of F is com- 
pared to the experimental ones [53| hi Fig. [T3] The 
theoretical and experimental curves show the same two- 
peak structure, but the calculated peaks appear at en- 
ergies which are about 3 MeV too low. This difference 
with the experimental dipole response is also found in 
the calculations by T. Inakura et al. (25[. In fact, the 
experimental mean energy that we extract from the data 
in the range [E min , £ max ] = [15, 30] MeV is equal to 22.1 
MeV, to be compared with our value of 19.8 MeV. The 
Gogny-force calculation by S. Peru et al. [H[ instead dis- 
plays the peaks at about the observed energies, with a 
mean energy of 23.0 MeV (cf. Table HI . 

One may comment that the excitation energies of the 
lowest peak predicted by our calculation (17 MeV) is just 
above the threshold for particle emission. At these ener- 
gies, a considerable part of photon absorption events will 



lead to a 7-ray cascade rather than particle emission. In 
fact, the 24 Mg(e, e') inelastic scattering experiment by 
Titze et al. [54[ displays a peak around 17 MeV (split 
into two rather narrow components), whereas the (7, n) 
and (7,p) data only display a tiny and statistically in- 
significant peak at this energy. On the other hand, (e, e') 
reactions are more complicated to analyze theoretically 
in a precise way than hadronic inelastic processes and 
photon absorption, since a consistent calculation of the 
cross section should go beyond the long wavelength ap- 
proximation, and also take into account magnetic inter- 
actions. 

For 26 Mg, the fractions of the EWSR may be com- 
pared to three experiments (55rl57l |. Fig. [TH displays the 
comparison between our curves and that of Fultz et al. 
[56j who, respect to the more recent work of Ishkhanov et 
al. [13], also measure the 26 Mg(7,pn) cross section. The 
three experiments agree with each other with respect to 
the overall width of the El-response, and they all dis- 
play two maxima at 18 MeV and 22 MeV separated by a 
shallow minimum around 20 MeV. This is in qualitative 
agreement with our calculated curves, especially when 
the THO basis is used and the QRPA spectra is folded 
with a Lorentian distribution with T=3 MeV. 

Combining the information from Table II and III, one 
notices that calculations based on a given interaction pre- 
dict a rather small shift of the mean energy of the giant 
dipole resonance, when comparing the two nuclei, 24 Mg 
and 26 Mg. This is at variance with experiment, accord- 
ing to which the mean energy shifts down by 1.7 MeV. 
In this way, a calculation which is in accordance with 
the data for 24 Mg, will disagree with data for 26 Mg, and 
viceversa. 

For completeness Fig. [15] shows the effect of leaving 
out the spin-orbit and Coulomb part of the interaction 
and one sees that this has only a minor effect contrary to 
the considerable effect obtained for the quadrupole mode. 
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FIG. 13. Fractions of IV dipole EWSR in 24 Mg, in HO (upper 
panel), and in THO (lower panel) basis. Dashed and dot-dashed 
curves are obtained from folding of QRPA spectra with a Lorentian 
distribution having T=l MeV and T=3 MeV respectively. These 
are compared with the experimental data by Irgashev et al. (solid 
curve) . 
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FIG. 14. Fractions of IV dipole EWSR in 26 Mg, in HO (upper 
panel), and in THO (lower panel) basis. Dashed and dot-dashed 
curves are obtained from folding of QRPA spectra with a Lorentian 
distribution having T=l MeV and T=3 MeV respectively. These 
are compared with the experimental data by Fultz et al. (solid 
curve) . 



B.3. IS + modes 

In Fig. [16] we display the experimental fraction of 
EWSR of Youngblood et al. [HJ and those obtained in 
the present work for the IS monopole mode of 24 Mg. The 
theoretical response shows a two-peak structure, with a 
low energy peak well defined around 18 MeV and a frag- 
mented high energy component which appears to be more 
sensitive to details and to the choice of the basis. In- 
deed, in the HO basis there is a significant contribution 
around 22 MeV, while in the THO basis this is in the 
energy range 25-30 MeV. Such peaks are not observed in 
the experimental curve, which is quite fiat and covers a 
broad energy interval. Except for the peaks, this spread- 
out behavior of the strength function is also found in the 
calculations. Also for the monopole mode, the calculated 



FIG. 15. The same as Fig. HOI but for IV 1 strength functions.. 



EWSR of 90 % is larger than the experimental value, 72 
± 10 %. From the similarity of the shape of the curves 
one can expect that the calculated centroids and the ex- 
perimental ones are in good agreement, as given by the 
numbers of Table [H] 

Fig. IT7l shows the calculated fraction of EWSR for the 
IS monopole mode of 26 Mg. As for 24 Mg, the behavior 
of the response is quite sensitive to the choice of the ba- 
sis. The curve displayed in the HO basis has a behavior 
similar to that obtained by S. Peru et al. i.e. the 
resonance has a prominent peak around 20 MeV and a 
broad component at higher energy. In the THO basis 
the low-energy peak of the resonance is shifted down by 
1.5 MeV and there is a significant contribution around 
26 MeV. 
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FIG. 16. The same as Fig. [9] but for fractions of IS monopole 
EWSR. 



Fig. [TS] shows the calculated IS quadrupole strength 
functions for the neutron-rich nucleus 34 Mg in the HO 
and THO basis. One can see a low-lying peak at 
2-3 MeV and a giant resonance at 15-22 MeV. The 
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FIG. 17. Fractions of IS monopole EWSR in 26 Mg, in HO basis 
(dashed curve) and THO (solid curve) basis. 



FIG. 18. IS 2+ strength functions in 34 Mg, in the HO (upper 
panel) and in the THO (lower panel) basis for the Q n =Q+ (solid 
curve), 1 + (dashed curve) and 2+ (dashed-dotted curve) excita- 
tions. 



first low-lying state belongs to the £l v =0 + component 
and it is mainly constructed by the neutron pp excita- 



tions [N, n 
3,3,0,±± 



2,0,2,±f + 



A,±fr] = 

2 

, the first of them coming from 

2 



3,2,1,±§ + 



(W 3/2 ) S 



the latter two from (I/7/2J • The features of these tran- 
sitions are shown in Tables [TV] E] for the HO and THO 
basis respectively. These results are in fair agreement 
with those of Table I of ref. (36J. One may notice that 
the two approaches give basically the same results and 
that the difference of about 200 keV between the two en- 
ergy peaks (2.64 MeV in the HO basis and at 2.48 MeV 
in the THO basis) is due to the unperturbed 2qp transi- 
tions with energy E qp x + E qp x . From the values of these 
2qp energy transitions one may conclude that the main 
contribution of the residual interaction is to shift down in 
energy the IS quadrupole response, rather independently 
on the choice of basis. 

In Fig. [19]we show the total isoscalar 2 + response func- 
tions in the HO and THO basis. When the spin-orbit and 
Coulomb residual interaction are omitted, the low-lying 
peak belonging to the n w =0 + , 2 + components is shifted 
up in energy by about 0.2 MeV. The giant resonance is 
instead shifted down in energy by about 0.6 MeV respect 
to the value obtained in the fully consistent calculation. 
The behavior of the strength functions calculated with- 
out spin-orbit and Coulomb interactions is quite close to 
the corresponding calculations by K. Yoshida et al. [35j ] - 
The giant resonance of our calculations is instead shifted 
up in energy by about 1.3 MeV respect to that in ref. 
[35(. This is coherent with the former comparisons be- 
tween our ISGQR and those of K. Yoshida et al. in 20 O 
where a difference of about 0.6 MeV is found, and in 24 Mg 
where there is a difference of 1.8 MeV. It seems that the 
blue shift displayed by the peaks in Fig. [19] as compared 
to those obtained in the calculations of K. Yoshida et 
al. without taking into account spin-orbit and Coulomb 
contributions in the residual interaction, increases with 
the intrinsic deformation of the nuclear system. 
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FIG. 19. The same as Fig. [TO] but for 34 Mg. 



V. CONCLUSION 

We described the development and testing of a fully 
consistent QRPA method to calculate linear response of 
axially-symmetric-deformed nuclei employing the canon- 
ical HO or THO HFB basis. The same Skyrme force is 
used in both HFB and QRPA approaches in all ph, pp 
and hh channels. The method is applied to study the 
responses of 20 O, 24_26 Mg and 34 Mg. 

For 20 O we showed the effective role of self-consistency 
in the IS 2 + and IV 1~ responses. Within this context, 
we performed calculations with and without the spin- 
orbit and Coulomb terms of the residual interaction. One 
can conclude that for this spherical nucleus these terms 
mainly act on the giant resonances, shifting their centroid 
up in energy by several hundred keV. 

We carried out detailed studies of the IS quadrupole, 
monopole and IV dipolc responses of deformed prolate 
and oblate 24 ~ 26 Mg. A microscopic analysis of the low- 
lying IS 2 + vibrations showed, for the open shell nucleus 
26 Mg major contributions to the response function aris- 
ing from pp transitions, while in the case of the nucleus 
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KK [2,0,2,±| + ] [3, 2, 1, ±§ J' [3,3,0,±i ] 





n 


n 


n 


v\ 


0.25 


0.68 


0.92 


£K 


-3.18 


-4.91 


-6.84 


Xkk 


-0.64 


0.65 


0.29 


Ykk 


-0.08 


0.03 


0.04 


E qp K + EqpK 


3.79 


3.72 


5.96 



TABLE IV. QRPA amplitudes [N, n z , A, ±0?]%- with isospin 
q (q = n for neutrons, q — p for protons) of the IS quadrupole 
Q. 7T = 0+ mode at 2.64 MeV for the HO basis in 34 Mg. The 
quasi particle occupations v%, the single particle energies ek 
and the 2qp excitation energies E qp k + E qv K , are given. Only 
components with X K k ~ Y K k > 0.01 are listed. 



KK [2,0,2,±| + ] [3,2, 1,±| ] [3,3,0,±i ] 



<IK 


n 


n 


n 


v\ 


0.24 


0.70 


0.92 


SK 


-3.20 


-4.91 


-6.84 


Xkk 


-0.64 


0.67 


0.27 


Ykk 


-0.09 


0.03 


0.04 


EqpK + EqpK 


3.58 


3.49 


5.72 



TABLE V. The same as Table [TV] but for the THO basis at 
2.48 MeV. 

24 Mg only ph (RPA) states are found to be important. In 
the energy region of giant resonances, the role of deforma- 
tion manifests itself in the splitting between the different 
projections which is more pronounced for the strongly de- 
formed 24 Mg (inhomogeneous damping). On top of this 
splitting comes a fragmentation, especially of the high- 
est energy projections of the giant modes. These effects 
lead to broad strength distributions. In this respect, our 
calculations confirm earlier theoretical results, and our 
consistent inclusion of the spin-orbit and Coulomb parts 
of the residual interaction only introduces minor changes 
concerning the giant modes, shifting the centroid by typ- 
ically 1 MeV and leaving the widths rather unaffected. 
The calculated strength functions in regions of giant vi- 
brations are compared to the available data, yielding gen- 
erally an overall account of the experimental findings. 

For 34 Mg, we give a description of the total IS 2 + re- 
sponse concentrating on the microscopic structure of the 
low-lying states which are in overall accord with the the- 
oretical results of K. Yoshida et al. [H, Hfl- 

We plan to optimize the present code, to be able to 
perform systematic calculations for both spherical and 
deformed systems in particular on the isotopes of Mg and 
of O. The aim is to extend our analysis to light drip line 
nuclei, so as to be able to study also pymgy resonances 
as well as other collective states typical of these exotic 
species. From a technical point of view, we also plan to 



assess for these exotic nuclear systems the significance 
of applying the THO-basis, which should more properly 
take into account the extended tails of wave functions 
than the HO basis. 

In a general perspective, the present work represents 
the first, unavoidable step for a consistent and more sys- 
tematic study of collective modes in nuclei, in particu- 
lar light exotic nuclei, taking properly into account also 
medium polarization effects. A study of the core po- 
larization effects in Al isotopes has been performed by 
K. Yoshida [38} by employing the quasiparticle- vibration- 
coupling model on top of the deformed HFB plus QRPA 
using a Skyrme interaction. 

The plan is to use the resulting states to calculate the 
role of the exchange of phonons between nucleons mov- 
ing in time reversal states close to the Fermi energy has 
in Cooper pair binding in exotic, deformed (as well as 
spherical) nuclei. 
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APPENDIX: EVALUATION OF INTERACTION 
MATRIX ELEMENTS 

A. Canonical wave-functions 

In this appendix we discuss the procedures used to 
evaluate some of the characteristic ph effective Skyrme 
matrix elements for the canonical wave functions having 
the general form: 

<$> K (r,cr) 




(10) 

with spin-up and spin-down components (pKt{ r -i-i z )i 
PKi{r_L, z) obtained in the present work in the harmonic 
oscillator basis, i.e. by the associated Laguerre and Hcr- 
mite polynomials [58j . The general form of these matrix 
elements and their effective Skyrme parametrization is 
given by J. Terasaki et al. Q in Eqs. (B12) to (B19). 
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Also, for the QRPA formalism we refer to J. Terasaki et 
al, Eqs. (A1)-(A6) in ref. QJ]. 

Concerning the most simple interactions, namely those 
proportional to the contact function 6(1 — r'), one has to 
evaluate the overlaps, 



er • er' = — \a+a' 

2 L + " 



(13) 



The algebra containing the 6(r— r 1 ) function is the same 
as above, leading to the result 



(KK'\5{r- t*)\LL') 



= j dr [$* K (r,tr)$ L (r,tr)] \®K' ( r i <T )*i' ( r > <T )] 



(11) 



overlaps which in terms of the canonical wave-functions 
(□H read 



(ifi^Xr- r')|Li'> - — 6n K+ n K „n L+ n L , 



x / r ± dr ± dz [ip Kt (r ±1 z)ip Lt (r ±1 z) 
+¥>Kl(r±, z)cp L i(r±, z)] 

x [(pK"t( r -L> z ) ( PL't(r±,z) + ipK'i(rj_,z)tp L >i(rx, z)] . 

(12) 

Next, we turn to the interaction er • cr'5(r— r J ). Here, 
the spin product is written in terms of lowering and rais- 
ing spin operators er+ and er_ 



(KK'\a ■ a'5{r-r')\LL') = — <Sn Jf +n A -„n i +n I 



x / r ± dr ± dz {[ip Kt (rj_, z)ip Lt (r±, z) 
- ( PKi(r±,z)ip Ll (r ± ,z)} 

x [<PK>f(r±, z)(p L , t (r ± , z) - (p K , l (r ± ,z)tp L , l (r ± ,z)} 
+ 2 [¥>jrt(rL, ^)^if'^(r_L, z)ip L i{r±, z)ip L , t (r ± , z) 
+ ^A'i(r_L, z)cp K '-t(r_L, z)(p Lt (rj_, z)tfL'i(r±, z) \ , 

where one can recognize the action of <J z a' z in the first 
term, and likewise h<?+cr'_ and ^o_a' + in the two subse- 
quent terms. 

It turns out that the momentum dependent parts of the 
residual interaction, involving differentiation in terms of 
the k and operators, give rise to quite many terms. 
Expressing one of these in full length 



k 2 = -iv 2 - -V' 2 + -V-V, 
4 4 2 



(15) 



one obtains 
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(KK'\8(r- r J )k 2 \LL') 



= -^u K +a K ,,n L +n L , / r ± dr ± dz 



<PKf{r±,z) 



1 d_ 

r± dr± dr 2 . 



(n 



3 + ^\m(.r±,z) 



1 d d 2 
r± dr± dr 2 . 



<9z 2 



x [ipK'-\-{rx,z)<PL'T-(rx,z) + ip K 'i(r_L, z)ip Ll i(rj_, z)] 

- -^Sa K+ Q K , tnij+ Q L/ J r^dr^dz [<p K ^(r±, z)ip L ^(rx, z) + <PKi(r±,z)<p L i(r±,z)] 



<PK't(r±,z) 



+ <fK'i(r±,z) 



1 d 



d 2 



fa - 1: 



rj_ dr± dr 2 . 



d 2 \ 

-^Avl^uz) 



1 d 



d 2 



fn, 



dz 2 



1 



r ± dr±dz I ip Kt (r ± , z)ip K ^(r±, z) 



¥>L>i{r±,z) 



(16) 



dr i 



+ (fiKt( r ±, z ) l PK'i(r±,z) 



'd<p L ^r ± ,z) 9 m (n,z) (fix - |) (£V + |) 



dri 



dr I 



9rj 



+ 



dtpL^(r±,z) dy£'i(r_L,z) 

fa + h) fa-h) 



+ ifiKi{r±,z)ip K 't(r±, z ) 



dipLj.(r±,z) d(p L ^(r±,z) 



¥>Ll(r±, z)ip L > t (r±, z) + 



dr± dr^ 
dfLi(r±, z ) d(p L > t {r ± ,z) 



+ <PKi.(r±,z)<p K >i(ri.,z) 



dz dz 
dtpLi{r±,z) d<p L >i(r±,z) (fix + |) (fti< + §) 



9rj 



dr± 



VLi{r±, z )(pL'i(r±, z ) 



+ 



difLi{r_L,z) dip L >i(r ± ,z) 



dz 



dz 



B. Spin-orbit interaction 



I 

parts, as they can be expressed in terms of raising and 
lowering operators as 



i<j> 



d 

dr± 



and 



The spin-orbit term i(cr + cr') ■ k? x S(r— r J )k of the 
residual interaction is likely the most involved to eval- 
uate. However, it becomes conceptually simple when 
it is interpreted as a volume product of three vectors 
A ■ (B x C). Next, the cartesian components of the vec- 
tors are replaced by the components of the spherical ten- 
sors of rank 1 (59|. This bears some resemblance to the 
treatment in the spherical case [l4| • In the deformed case, 
the spin spherical tensors raising and lowering operators, 
and the spherical tensors of the differential operators also 
acquire a more simple form than their cartesian counter- Still, the volume product contains six terms, namely 

I 



.on 



j_d_ 

rj_ dcf. 



i d 
r± d(j) 



(17) 



(18) 



(KK'\i{(7 + a') • Jtf x 8{r- r )fc|£L'^ = i [{KK'\ (o+ + a' + ) (V z - V'J S(r- r ) (V_ - V'_) \LL') 

- (KK'\ (a++a' + ) (V_ - V_) S(r- r') (V z - V'J \LL') + (KK'\ (a z +a' z ) (V- - V'_) S(r-r') (V+ - V'+) \LL') 

- (KK'\ (a z + o' z ) (V+ - V' + ) S(r-r') (V- - V'_) \LL') + (KK'\ (a. + a>_) (V+ - V^) 5(r- r) (V, - Vj \LL') 
(KK'\ {a. + a'_) (V z V J 5{r r') (V+ - V^) \LL')] . 

(19) 
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In the above relation, the operators to the left of 6(1 — r 1 ) 
act on K, K', while the operators to the right of 6(1 — r 1 ) 
act on L, L' . The sum over the six terms of Eq. (fT9|) 
corresponds to a sum over the six permutations of the 
set {+, z, -}. 



C. Coulomb interaction 

Due to the short range character of the nuclear residual 
interaction, which is expressed through the S(r— r 1 ) func- 
tions, the four-dimensional integrals J r±dr±dzr' ± dr' ± dz' 
are directly replaced by two-dimensional integrals. How- 
ever, for the direct term of the Coulomb interaction 

^d/rect = \r-r J \ one neeo ^ s to carry out the full integra- 
tion. In order to exploit the cylindrical symmetry, we 
make use of an expansion method applied to astrophysi- 
cal problems by Cohl and Tohline [Ho, EH 



r'l 



'\fr± 



OO 

■ £ 

rn— — 00 



Q m 1 (x-)e l 



i(4>-4>') 



(20) 



(KK'\ V d e l f ect \LL') = 2e 2 (2tt) 2 S m ,n K -n L • 6 m ,n L ,-n K , 

00 OO OO CO 

x J r^dr^ J r'^dr'^ j " dz J dz' 

00 

x [<PKf(rx, z)v>Lt(rx, z) + VKi(r±, z)ip Ll (r ± , z)} 

Q ro _i (X-) + KK> ■ 7TL' • (-I)™ Q m -i 



iry/r±r' ± 

x [^'t( r l' ^)fL't(r'±, z') + <PK>±(r±,z , )<PL>±(r'x,z')] , 

(24) 



with 



X+ 



r\ + r'l + {z + z') z 



2r 1 r', 



(25) 



Since itk> -ttl' ■ (— l) m does not depend on the projections 
Q.K — f^L i ~ > this factor is the same for spin-up 
and spin-down wave functions. 



where 



Here Q m _L is a Legendre function of the second kind 
of half- integer degree [58[ . We have checked the conver- 
gence of this expansion, which is quite rapid when the 
parameter \ is not too close to 1. Inserting the canonical 
wave- functions, the integrals over the azimuthal angle se- 
lect the order m of the projection of the angular momen- 
tum on the z-axis, such that m = fix — = fir,' — Qk 1 ■ 
The functions Q m _i can readily be evaluated by sim- 
ple integrals and tabulated over a suitable range. In the 
calculation of the matrix element 



(KK'\vg* et \LL') = JJdrdS 



x $^(r, tr)$* K , {r / , a)- -77 $l(t", <t)$ l , {r> , tr), 

\r—r\ 

(22) 

it is important to take into account the symmetry of the 
wave- functions 



$K(-r f tr) = ir K ®K(r,er), (23) 

7tK — ±1 being the parity of the state K depending on 
the sign of z. By exploiting these symmetries one needs 
only to carry out integration over z,z' > 0, aside from 
integration over the two angles <fi and <f>' 
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